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Giant Flare of SGR 1806-20 from a Relativistic Jet 
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Abstract 

Japanese magnetospheric explorer GEOTAIL recorded a detailed light curve during the initial 600 msec 
of a giant flare from SGR 1806-20 on December 27, 2004. We show that the observed light curve is well 
explained by an emission from relativistically expanding fireballs, like those of gamma-ray bursts (GRBs). 
Especially, the observed rapid fading after 500 msec suggests that eject a is collimated in a jet. We derive 
an upper limit on the jet opening half-angle of 0.2 radian that is as narrow as those of GRBs. 
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1. Introduction 

Soft gamma-ray repeaters (SGRs) are most likely highly 
magnetized neutron stars, so called magnetars (Thompson 
& Duncan 1995). On December 27, 2004, a giant 
flare from SGR 1806-20 illuminated the Earth (Boggs et 
al. 2004; Hurley et al. 2005; Palmer et al. 2005) with 
the gamma-ray flux more than ~ 10 6 times that of the 
typical cosmological gamma-ray bursts (GRBs) that are 
the most violent explosions in the universe. The giant 
flare has an initial spike lasting about 600 msec with an 
isotropic equivalent energy of ~ iq 46-47 ergs, that is fol- 
lowed by a pulsating tail lasting 400 seconds with energy 
of ~ 10 44 ergs. Only two giant flares from SGRs had been 
recorded before this flare: they occurred on March 5, 1979 
from SGR 0526-66 (Cline et al. 1980) and August 27, 1998 
from SGR 1900-14 (Hurley et al. 1999), respectively. The 
initial spike of the most recent flare is ~ 10 2 times more 
energetic than the previous two events, while pulsating 
tails of the three events have comparable energies (e.g., 
Woods 2004). 

Because of the observed high flux density, most 7-ray 
detectors were saturated except for particle detectors such 
as GEOTAIL that successfully recorded a burst light curve 
in the brightest initial spike (Terasawa et al. 2005; see 
also Mazets et al. 2005). The burst was so bright that 
the light curve was clearly recorded down to three orders 
of magnitude below the peak flux. In the early epoch 
(t < 160 msec), the light curve is variable and mainly con- 
sists of two pulses. A gradual power-law like decay be- 
gins after the second peak (t > 160 msec) with a bump 
at t ~ 430 msec. This is followed by a rapid fading after 
t ~ 500 msec. Such a detailed light curve of the initial 
spike has been measured for the first time. It brings us 
a new key to understanding of the mysterious, poorly un- 



derstood SGRs. 

The detection of radio afterglows after giant flares 
suggests that SGRs eject relativistic outflows (Cameron 
et al. 2005; Gaensler et al. 2005; Frail, Kulkarni & 
Bloom 1999). A relativistic motion is also implied by 
the nonthermal flare spectrum observed by Mazets et 
al. (2005) and Palmer et al. (2005), otherwise pair for- 
mation occurs in a compact emission region, which makes 
thermal spectrum (Huang, Dai & Lu 1998; Thompson & 
Duncan 2001; Nakar, Piran & Sari 2005). Even if the spec- 
trum is thermal (Hurley et al. 2005), its hyper-Eddington 
luminosity implies a relativistic motion (Wang et al. 2005). 
The initial spike and the pulsating tail have different spec- 
tral features and temporal pulse profiles, which suggests 
that they have different origins: the initial spike and 
the radio afterglow may arise from the relativistic out- 
flow, and the pulsating tail may come from evaporating 
trapped fireballs. The kinetic energy of the outflow in- 
ferred from the radio afterglow is ~ 10 44-45 ergs (Cameron 
et al. 2005; Gaensler et al. 2005), which is much smaller 
than the isotropic equivalent energy of the initial 7-ray 
spike. Hence a collimated relativistic outflow is implied. 
In addition the isotropic equivalent energy of the initial 
spike is comparable to that of the exterior magnetic field 
£>~10 15 Gof the magnetar. Then the magnet ar cannot 
produce giant flares repeatedly (e.g., ~ 100 times) during 
its active time (~ 10 4 yrs) unless the emission is collimated 
or there is other energy reservoir. Therefore we have fair 
motivations to consider an anisotropic giant flare. 

If the giant flares of SGRs arise from relativistic col- 
limated outflows, they are similar to canonical GRBs 
from the cosmological distance (typically tens of Gpc). 
A sub-group of GRBs with long duration is thought to 
be caused by relativistic jets that originate in the col- 
lapse of a massive star (Hjorth et al. 2003; Stanek et 
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al. 2003). Energy is carried away from a compact source 
as kinetic energy of jets (Piran 1999; Piran 2004; Zhang 
& Meszaros 2004). This is converted into radiation by in- 
ternal shocks between shells, which make observed highly 
variable gamma-ray light curves, called prompt emissions 
of GRBs. Subsequently, at larger radii, outflow inter- 
acts with ambient circumstellar matter, producing exter- 
nal shocks, which are responsible for afterglows on much 
longer time scales in various wave lengths, such as radio, 
optical, and X-ray bands. 

In this Letter, we show that the observed light curve of 
the initial 600 msec is well explained by the emission from 
relativistically expanding fireballs, like those of GRBs. 
Especially, the observed rapid fading after 500 msec sug- 
gests that ejecta is collimated in a jet. We derive a robust 
upper limit on the jet opening half- angle of 0.2 radian that 
is as narrow as those of GRBs. 

2. Light curve of collimated outflow 

The light curve of the initial spike of the giant flare from 
SGR 1806-20 is very similar to the behavior in prompt 
GRB emissions (see Figure 1). We can interpret two 
pulses in the early epoch (t < 160 msec) as two inter- 
nal shocks (see below). The following decay is basically 
determined by the relativist ic kinematics, which is inde- 
pendent of the emission mechanism. Suppose a relativistic 
shell shines for a short period. Since the shell has a curva- 
ture, photons far from the line of sight (LOS) come later. 
The shell at higher latitude from the LOS has a lower ve- 
locity toward the observer, so that the emission becomes 
dimmer and softer as time goes because of the relativistic 
Doppler effect, which explains the observed power-law like 
decay during between 200 and 400 msec very well. If the 
emission is spherical (isotropic), however, such a decay 
should continue beyond t > 600 msec. This is inconsistent 
with the observation that the light curve rapidly fades af- 
ter t ~ 600 msec, which implies that the emission does 
not occur at larger angle from the LOS. In other words, 
the giant flare arises from a relativistic jet with a finite 
opening angle. 

In order to see the above arguments quantitatively, we 
consider a simple model for emission from a relativisti- 
cally moving jet which radiates photons when the shell 
is located at radius from ro to r e (Yamazaki, Ioka & 
Nakamura 2003). We assume that the cooling timescale 
is much shorter than other timescales, and hence consider 
an instantaneous emission at the shock front. We use a 
spherical coordinate system (t,r,0,(/)) in the Lab frame, 
where the 6 = axis points toward the detector at r = D, 
and the magnet ar is located at r = 0. The jet has an 
opening half-angle AO and a viewing angle V the axis of 
the emission cone makes with = axis. The emitting 
shock front moves radially from t = to and r = ro with the 
Lorentz factor 7 = 1/ y/l — f3 2 . Then the observed flux per 
unit frequency of a single pulse at the observed time T is 
given by 
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where f(v') represents the spectral shape, and 1 — 
P cos 6{T) = (c(3/r )(T - T ), 1 - (3 cos 0(t) = [1 - 
f3cos6(T)}/[(cp/r )(t-T )} and T = t - r /c(3. The value 
of A(j)(t) is 7T in the case of V < AO and < 0(t) <A6-0 V , 
while for the other case, it is given by 
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The normalization of emissivity A(t) is determined by the 
hydrodynamics. Here for simplicity we adopt the follow- 
ing functional form, 



A(t)=A 
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where the emission ends at t = t e and the released energy 
at each distance r is constant. Shell emits at radius from 
ro to r e = ftro, where K = t e jt§ > 1. The quantity Aq is a 
normalization constant. A pulse-starting time and ending 
time are given as 



Tgtart = T + (r /c/3)(l- 
T end = T + [(r /c/3)+t e 



/3cos(max[OA-A0])) , (4) 
-to](l-0CQB(6 v + AO)) ,(5) 

respectively. We adopt the following form of the 
comoving-frame energy spectrum, 

/(*/) = (^K) 1+Q exp[-(*/K)] , (6) 

where a is a power law index. Mazets et al. (2005) derived 
a nonthermal spectrum with a = —0.7 and an exponen- 
tial cut-off at 800 keV, while Palmer et al. (2005) gave 
a nonthermal fit with a = — 0.2 and a cut-off at 480 keV. 
Although Hurley et al. (2005) reported a black-body spec- 
trum (a = 1) with a temperature higher than ~ 240 keV, 
at least a portion of the giant flare may be nonthermal 
because the spectrum was determined with low time res- 
olution. A previous giant flare in SGR 0526-66 may also 
have had a nonthermal spectrum (Fenimore, Klebesadel & 
Laros 1996). In addition, some giant flares from SGRs in 
nearby galaxies (within 40 Mpc) may have been detected 
by BATSE as short GRBs whose spectra are likely non- 
thermal (Palmer et al. 2005). In this Letter, we assume 
a nonthermal spectrum, which is naturally produced by 
shocks in a relativistic outflow. 

Equations (1), (3) and (6) are the basic equations to 
calculate the flux of a single pulse, which depends on fol- 
lowing parameters: j0 v , jAO, jv' , ro/c/3j 2 , a, to, t ei and 
r 2 j 2 A /D 2 p. We fix j6 v = 0, jhv' = 800 keV, and a = 
—0.6 in the following. Then we find out best fit values for 
other parameters: 7Ar9, ro/c/?7 2 , and n (=t e /to = r e /ro). 
We find that these parameters do not depend on a and 
7^ so much, even if we assume a black-body spectrum. 

Figure 1 shows the result. The fit is surprisingly good 
considering the very simple model. The second pulse, 
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which has a duration of Tab ~ 50 msec, and the associ- 
ated power-law like decay lasting Tbc ~ 500 msec, are well 
fitted by our model with jAO = 3.0, ro/7 2 = 2.6 x 10 8 cm, 
and n = 12.5. We have checked that the uncertainty com- 
ing from the spectral shape is at most a factor of 2. We 
also show an example of the theoretical modeling for the 
first pulse (t < 80 msec) by black dotted line in the same 
figure. In this case the opening half-angle is not well con- 
strained because the power-law like decay is masked by 
that of the second pulse. We can make a similar modeling 
for a bump at t ~ 430 msec, though it is not shown in the 
figure. 

The opening angle of the jet is constrained by the light 
curve in a kinematical fashion (see Figures 1, 2). The 
duration of the brightest epoch (80 < t < 130 msec) is 
determined by the crossing time of the shell through the 
emitting region (ro < r < r e ) as 

(r e -r )(l-/3) 



Tab = 



50 msec. (7) 



c(3 v " ^2c 7 2 

On the other hand, the following power-law like decay 
lasts for ~ 500 msec, which is approximately given by the 
angular spreading time 

r e (l -cosA6>) 



Tbc ■ 



2c 7 2 



500 msec. (8) 



In other words the wider the jet, the later the onset of the 
steep decay. Eliminating ro/2cj 2 from these equations, 
we derive (jA6) 2 ~ 10(1 - ft" 1 ) < 10 and hence jAO < 3. 
The uncertainty is at most a factor of 2. Furthermore, 
combining with 7 > 25 required to avoid pair formation 
(Nakar, Piran & Sari 2005), we obtain a firm upper limit, 
AO ^0.2 radian, which is very similar to those of GRB 
jets inferred from the observed break in afterglow light 
curves (Harrison et al. 1999). In particular, we disfavor 
models with isotropic emission for the beginning epoch of 
the giant flare. 

3. Discussions 

We have shown that the initial spike of the giant flare 
SGR 1806-20 is well explained by emission from a rela- 
tivistic jet directed toward us, and that the opening half- 
angle of the jet is less than 0.2 rad. Since the isotropic 
equivalent energy of the giant flare is 5 x 10 46 ergs with 
assumed distance of 15 kpc (Terasawa et al. 2005), the 
collimation-corrected energy is less than 5 x 10 44 ergs for 
AO < 0.2 and the flare is rather economical than previously 
thought. This may alleviate an extreme situation that an 
isotropic flare demands almost all energy of dipole mag- 
netic fields of SGR 1806-20. The size and the light curve 
of the radio afterglow from SGR 1806-20 also favor smaller 
energy ~ 10 44 -10 45 ergs (Cameron et al. 2005; Gaensler et 
al. 2005) than the isotropic equivalent energy of the flare, 
which also suggests a jet opening half-angle AO ~ 0.2. 

External shock scenario is unlikely because it is difficult 
to keep high Lorentz factor and have a steep decay at 
500 msec. The third bump at t ~ 430 msec may be caused 
by an additional internal shock, while other reasons such 
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Fig. 1. Comparison of theoretically predicted light curves 
with the observed data. Green and blue points are 
background-subtracted MCP and CEM data, respectively 
(Terasawa et al. 2005). The second pulse and the follow- 
ing power-law like decay are modeled by red lines, which 
have AO = 27 -1 ,37 -1 , and 47 -1 from left to right with 
ro = 2.6 x 10 8 7 2 cm and r e = 12.5ro- The rapid fading at 
t ~ 600 msec is most consistent with an opening half-angle of 
A6 = 37" 1 . The black dotted line shows an example of the 
theoretical modeling for the first pulse (t < 80 msec). 
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Fig. 2. A schematic picture of the jet emission. An observer 
resides far on the right side. A thin shell emits gamma-rays 
while it crosses the hatched region (ro < r < r e ). Each red 
arrow represents the emitted photon at each place. The ob- 
served duration ~ 50 msec of the second pulse (80 < t < 130 
msec) in Figure 1 is determined by the shell crossing time, 
Tab, i-e., the difference of the arrival time of two photons A 
and B that are emitted when the shell crosses radii ro and 
r e , respectively. The observed duration of the power-law like 
decay after the second pulse (130 < t < 600 msec) in Figure 1 
is determined by the angular spreading time, Tbc, i- e -> the 
difference of the arrival time of two photons B and C emit- 
ted simultaneously. The wider the jet is, the later the rapid 
fading begins. 

as inhomogeneities on the jet are possible since its peak 
flux is much smaller than the others. 

A relativistic jet begins sideway expansion when the jet 
Lorentz factor becomes 7 ~ (A#) -1 . This epoch is ob- 
served at Tj e t ~ 9 min 

(£45/S)) 1/3 (A<9/0.1 r ad) 2 , where 
£?45 and no is the total energy confined in the jet and 
the number density of the ambient matter. After that 
time, the jet decelerates abruptly. Hence in the epoch 
of radio observations (6-20 days after the flare), the 
outflow became Newtonian or sub-relativistic and nearly 
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isotropic, which is consistent with the radio observation. 
Nevertheless some degree of anisotropy may remain and 
produce the observed elliptical image and polarization of 
the radio afterglow (Gaensler et al. 2005). 

It has been discussed that giant flares from other 
SGRs resemble classical GRBs in spectroscopic charac- 
ters (Fenimore, Klebesadel & Laros 1996). This possibil- 
ity is strengthened by our present result that the recent 
giant flare of SGR 1806-20 is a jetted emission like GRBs. 
However, in our scenario, there should be many more 
misaligned SGRs, which will show up only in isotropic 
emission. If the pulsating tails are isotropic emission, 
there should be many events consisting of only pulsating 
tails, though such events that are bright enough to trigger 
e.g. BATSE have not yet been reported. One possibility 
to resolve this "statistical" problem is that the pulsat- 
ing tail is also collimated. According to the magnetar 
model, pulsating tail emission arises from an evaporating 
trapped fireball (Thompson & Duncan 1995; Thompson 
& Duncan 2001). Its size can be comparable to the radius 
of the magnetar. However, as discussed in Thompson & 
Duncan (2001), due to the QED effect, the radiative trans- 
port across the magnetic field lines is concentrated to the 
foot-point at the magnetar surface, which may be a reason 
for the collimation of the pulsating tail. The collimation 
angle would largely depend on the magnetic field config- 
uration. However, the broad pulse profile (larger than 
1 sec) of the tail disfavors the narrowly collimated pulsat- 
ing tail emission. Another way to resolve the problem is 
to introduce some envelope around the main jet. In the 
GRB case, it is now widely argued whether the angular 
structure of the jet is uniform, Gaussian, power-law, or 
two-component (Zhang & Meszaros 2004). Similar to the 
GRB, it may be possible that the SGR jet also has a cen- 
tral core with AO ~ 0.1 and an envelope with wider solid 
angle that can produce less energetic flares. Indeed, the 
past two giant flares show much smaller isotropic equiva- 
lent energies, ~ 10 44 ergs, while all three flares had pulsat- 
ing tails with a similar energy ~ 10 44 ergs. If the central 
core is seen off-axis and the observer points to the enve- 
lope, such less energetic giant flares can be observed. At 
that time, the intense emission from the central core may 
be negligible because of the relativistic beaming effect. If 
such an envelope has a wide solid angle, the statistical 
problem can be resolved. Indeed, exponentially decaying 
tail after t ~ 600 msec recorded by B AT / Swift (Palmer et 
al. 2005) may be such an envelope emission. 

The relativistic jet may not be generally launched along 
the rotation axis because the magnetic energy domi- 
nates over the rotational energy. The duration for the 
magnetar to eject a jet is less than the duration of a 
pulse, tdur < Tab ~ 50 msec. Thus the rotation angle 
during the jet ejection is at most $ rot ~ 2tt (td U r/T p ) < 
0.04(^ nr /50 msec) rad, where T p = 7.56 sec is the ro- 
tational period of SGR 1806-20 (Hurley et al. 2005). 
Therefore, the overall opening half-angle is at most AO + 
d rot /2. Since this is comparable to A0, the effect of the 
rotation of the magnetar may not be so large but may 
still affect the jet structure and generate the envelope as 



discussed in the previous paragraph. 

If the energy source of the flare was a dipole outer mag- 
netic field and ~ 10% of its energy was converted to 7-rays 
as in the case of the isotropic flare, the spin down rate, 
P oc B 2 /P, where P is the rotational period of the mag- 
netar, would change by ~ 10% after the flare. Thus the 
isotropic flare might be tested by the observations of the 
spin down history. 

It is widely believed that anomalous X-ray pulsars 
(AXPs) are the same kinds of objects as SGRs (e.g., 
Gavriil et al. 2004). Kulkarni et al. (2003) proposed the 
main difference between AXPs and SGRs is the time- 
dependent geometry of the magnetic fields, i.e., AXPs are 
older and less active than SGRs. Here we suggest that the 
viewing angle may also contribute to the different appear- 
ances of AXPs and SGRs. If we see the passive (active) 
part of the magnetar surface, the magnetar would show 
soft (hard) persistent X-ray emissions and look like AXP 
(SGR). This might be consistent with no giant flares from 
AXPs because the energetic jet core responsible for the 
giant flare is always misaligned. 
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